Quasiprobability methods for multimode conditional optical gates. 
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We present a method for computing the action of conditional linear optical transformations, 
conditioned on photon counting, for arbitrary signal states. The method is based on the Q-function, 
a quasi probability distribution for anti normally ordered moments. We treat an arbitrary number 
of signal and ancilla modes. The ancilla modes are prepared in an arbitrary product number state. 
We construct the conditional, non unitary, signal transformations for an arbitrary photon number 
count on each of the ancilla modes. 
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It has been known for some time that non deterministic non linear transformations are possible with linear optical 
networks when some subset of the input modes (ancillas) are prepared in single photon states before the optical 
network and directed to photon counters at the network output [l|, 0, 0, [4, 5]. The conditional state of all non-ancilla 
modes (the signal modes), conditioned on a particular count on the output ancilla modes, is given by a non unitary 
transformation of the input signal state and can simulate a highly nonlinear optical process. This transformation is 
defined in terms of a conditional measurement operator acting on the signal modes alone. It is determined entirely 
by the total linear optical unitary transformation (acting on both signal and ancilla modes), the known ancilla 
input number states and the ancilla photon counts at the output. A number of approaches have been proposed for 
constructing the conditional state transformation. Scheel et al.@ give an algebraic approach, while Lapaire et al.Q 
show how nonlinear such transformations can be. In Q van Loock et al. derive a set of criteria to decide whether 
a given measurement operator can be implemented. In this paper we show how quasi probability methods (the Q 
function) can be used to construct the measurement operator for arbitrary ancilla number states and arbitrary ancilla 
number counts. The method leads naturally to a diagrammatic representation of the measurement operator. 

Consider the situation depicted in figure [T] in this device N + K modes pass through a optical device that 
conserves total photon number. We will call this a linear optical system. The K ancilla modes are prepared in photon 
number eigenstates. At the output, photon number measurements are made on the ancilla modes alone. We seek the 
conditional state for the remaining N modes given the ancilla photon number count. 
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FIG. 1: A conditional linear optical gate. 

The linear optical device performs a unitary transformation on all the input states. It may be described by a 
unitary transformation of the form[^, [^] 



U(H) = eiq>[-ia?Ha] 



(1) 
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where 



a 2 



a N 



(2) 
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and H is a hcrmitian matrix. This transformation leaves the total photon number invariant, 

U f (H)a^ -aU{H) = a f • a 
It also induces a linear unitary transformation on the vector a as 

U\H)aU{H) = S{H)a 



(3) 



(4) 



One should not confuse the unitary transformation, U(H), acting on states, with the induced unitary representation, 
S(H), acting on the mode operators. 

The conditional state of the signal modes, \ip'} is then determined by 
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E(n\rh)\tp} s 



(5) 



where the observed count is represented by the vector of values m, and the probability for this event is p(m). The 
measurement operator is 



E(n\rh) = anc(m\U (H)\n) ar . 



with 



\m)anc = |mi)jv+i ® \m 2 ) n +2 ® ■ ■ ■ ® \m k )N+K 
In what follows we will need to know how coherent states transform under U (H) . We define 

\a) = \ati) ® \a 2 ) ® . . . |aAr +A -) 

where 
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(7) 
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Then 



U{H)\a) = \a(H)) = {a^H)) ® |a 2 (ff)) . . . |ajv+jr(5)) 



(10) 



where <J(i7) = S^(H)d , reflecting the well known fact that coherent states do not become entangled through a linear 
optical transformation of the kind considered here. 

To prove this result, we work with the Bargmann coherent states [Toj. 
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which arc related to Glauber coherent states by \a) 

U{H)\\a) 



exp[-at.d/2]||<3). Then 

U(H)e sf - s \0) 
U(H)e s - s tf(H)\0) 



10) 



\\a(H)) 



|0) 



where a(H) = S^i^d?. Now as oft{H).a(H) = eft .a, we can write U(H)\a) 
We now define the function 



\d(H)). 



U H {a) = (a\U(H)\a) 



i\\U(H)\\< 



(12) 



which completely determines the operator U (H) (which acts in the full tensor product Hilbert space for both signal 
and ancilla modes, as we now show. 

Due to the over completeness of the coherent state basis, this diagonal matrix element suffices to completely 
determine the operator To see this we note that if any operator is put in normal order 



the Q-symbol is defined by 



As a | a) = a | a) we see that 



i = ^ar im (a t )"a™ 



A(a) = (a\A\a) 



(13) 



(14) 



(15) 



It is easy to see that the Q-symbol for U (H) is given by 

Uh{o) = exp[a t 5(i?)a] exp [— eft ■ a] (16) 

The Q-symbol, A{a) completely determines the corresponding operator, A, by the following method. 

• Expand A(a) as a power series in (a*) n a™ 

• Order all conjugates to the left and make the replacements {a*) n i— > (<xj)" and a™ i— ► a" 

This mapping between Q-symbols and operators can always be done efficiently with respect to the number of modes, 
and can thus be efficiently implemented in a symbolic computer package. 
We can also define the Q-symbol for the measurement operator in Eq.([5]) as 



£ H {d s \n,m) 



where 



o. s . 



where the subscript s refers to the signal modes alone. 
Noting that 



1 



s (a s |£^(n|m)|a s ) s 

\ U n ) 

d 



(17) 



(18) 
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(19) 
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and that the partial derivatives commute with U(H), it is easy to see that 

£ H (a s \n,m) = [m N+ i\ . . . m N+K \n N+1 \ . . . n N+K l]' 1/2 [(d* N+1 )" lN+1 . . . (d* N+K ) mN+K (20) 
x (d N+l ) nif+1 . . . (d N+K ) nN+K exp[dtS(iJ)d]] 5 _ =o x e- s l' a > 

where 



ajv+2 

\ &N+K J 



(21) 



and 

d 



d* n = 

d n = 



da* 
d 
da„ 



This expression enables us to construct the relevant Q-symbol for the conditional state transformation by expanding 
the function as a power series and reading off the relevant coefficient. This function thus plays the role of 

a Feynman propagator in ordinary quantum mechanics, ft describes the many ways in which a classical coherent 
amplitude can propagate through the linear optical network, keeping track of the manyfold coherent amplitudes for 
reflection and transmission. This captures in a formal way a heuristic that has often been used to construct conditional 
gates in the past 

As an introduction we consider two modes, 04,0,2, coupled with a beam splitter interaction, so that the unitary 
transformation of the modes is given by 

*-(JS) (22) 

then the transformation is paramctcrised by a single parameter, 9 and 

01(6*) = cos 6a\ — sin 60,2 (23) 

02 {&) = cos 6a,2 — sin 6a\ (24) 

Then we see that 

U H (ai,a 2 ) = exp[(|ai| 2 + |a 2 | 2 ) cos6» + {a* 2 ai - a\oL 2 ) sin6i]e~ (|Ql|2 - |a2|2) (25) 

We now assume that photons are counted on mode a 2 and calculate the conditional state for mode a\ for the four 
cases: n = 0,1 and m = 0, 1. 

The conditional state of mode ai is given by, 

|V (n ' m) )i = r l— E(n\mM 1 (26) 
\JP\m\n) 

where 

p(m\n) = i(V'|i? t (n|m)£;(n|m)|V')i (27) 
The Q-symbol for the operator E(n\m) is 

d^J \da~2) [ ex P[(l a i| 2 + l a 2| 2 )cos0+(^a 1 -aJa 2 )sin0]] Q2=o e-l Ql l 2 (28) 
The corresponding operators are 

£(0| ) = ( C0S9 -^ n ( t)n fl n = . e -ln{co S 9)a\ ai . 

n=0 U ' 

E{1\1) = cos6'i?(0|0)-sin 2 6'al^(0|0)ai 
E(1|0) = -4sin^(0|0) 
E(Q\l) = sin 6*^(010)01 
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Now consider a three mode example defined by the transformation 

(sil S12 si 3 \ / ai \ 
S21 s 22 s 23 «2 (29) 
S31 S32 s 33 J \a 3 J 

We will regard 03 as the ancilla mode, prepared in a single photon state \n)3 and we will condition off a count of 
"1" at this mode after the optical system. Thus we have a three mode example with n 3 = 1, 7713 = 1 and we wish to 
calculate the two mode conditional measurement operator £7(1 11) acting on modes a±, a 2 . 
The formalism gives 

E(l\l) = s 33 A + (si3<4 + s 23 al)A(s 31 ai + s 32 a 2 ) (30) 
where the Q-symbol for A is given by 

A = exp[— (1 - sn)|ai| 2 - (1 - s 22 )\a 2 \ 2 + s\ 2 a\a 2 + s 2 ictia 2 ] (31) 

Using the theorem [l2| 

(a 1 ,a 2 \e Xia ^ a2 e^ ai+ ^ a2 e X2aia ^\a 1 ,a 2 ) 
= exp [-|a 2 | 2 (l - e^-) - |«i| 2 (l - e" 1 - A^e" 2 ) + Ka 2 Ai + a 1 a 2 A 2 )e' 12 ] 

we see that 

_ e ^ia{a 2 e fJ.ia\a 1 +fj.2ala2 e X 2 aia\ (32) 

with 



S22 



(33) 
(34) 

(35) 

(36) 

The form of this Eq. ([30|) suggests an interpretation as a 'sum over histories' for the different ways a photon can go 
from input to detection. For each of the terms in Eq. (|30p . we can assign a diagram, see figur^2] 
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1st term 2nd term 3rd term 4th term 5th term 

FIG. 2: A set of diagrams corresponding to each of the terms in S33A + aJsi34s3iai+aJsi3^4s32a2 + a2 S2 3 y ^ s 3 iai + a 2 S2 3 y ^ ,s 32 a 2- 
The modes are labeled 1, 2, 3 from top to bottom. 



As an explicit example, consider the situation where all beam splitters are 50/50. In this case 

1 



2V2 

The 



y/2 + 1 V2- 1 V2 

y/2-1 V2 + 1 -V2 ) (37) 
V2 -V2 -2 



E(l\l) = ~A + \{a\ - 4)A( ai - oa) (38) 
If the input state on modes a 12 is |?/>) = |1, 1), the (un normalised) output state is given by 

E(l\m, 1) = -|([2 ; 0) + [0,2» - -±=\1, 1) (39) 
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The norm of this state gives the probability for this event, 

P(l|l) = 5/16 (40) 

We have presented a method based on quasi probability distribution for constructing conditional linear optical 
transformations. The explicit construction involves a power series expansion of the Q-symbol for the linear optical 
network that defines the process. This method lends itself to a diagrammatic interpretation. More importantly it 
lends itself to a computer based symbolic computation. We believe this method will assist in finding new conditional 
state transformations. 
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